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Abstract 

This paper is concerned with the asymptotic spreading of a Lotka-Volterra co- 
operative system. By using the theory of asymptotic spreading of nonautonomous 
equations, the asymptotic speeds of spreading of unknown functions formulated by a 
coupled system are estimated. Our results imply that the asymptotic spreading of one 
species can be significantly fastened by introducing a mutual species, which indicates 
the role of cooperation described by the coupled nonlinearities. 
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1 Introduction 



In this paper, we consider the propagation of the following diffusion system 



< 




d±Aui(t, x) + riui(t, x) [1 — u\(t, x) + 61U2 (t, x)\ 
d 2 Au 2 (t, x) + r 2 u 2 (t, x) [1 - u 2 (t, x) + b 2 u\(t, x)] 



(1.1) 



in which ui(t,x),u 2 (t,x) denote the densities of two collaborators at time t > and 
location x S R in population dynamics, all the parameters are positive and b\b 2 < 1 such 
that (jl.ip has four spatial homogeneous steady states (for short, four equilibria) 



It is well known that (k\, k 2 ) is asymptotic stable while (0, 0), (1, 0) , (0, 1) are unstable in 
the corresponding spatial homogeneous system of (jl.ip . 

Recently, Li et al. [11] have investigated the traveling wavefronts of (jl.ip by using 
the theory established by Weinberger et al. [27J, and the authors proved that the mini- 
mal wave speed of traveling wavefronts of (jl.ip can be linearly determinate (see [3lll9|). 
In the modeling of population invasions (see Shigesada and Kawasaki [23] for many im- 
portant historic records), the linear determinacy indicates that the minimal wave speed 
can be formulated by the parameters appearing in the system linearized at the invadable 
equilibrium which often is unstable in the corresponding kinetic system. In population 
dynamics, besides the minimal wave speeds of traveling wavefronts, the asymptotic speeds 
of spreading may also be linearly determinate, especially for the scalar equations, we refer 
to Aronson and Weinberger [I], van den Bosch [3], Diekmann [TIE] 5 Hsu and Zhao |10| . 
Lui |16jll7j. Mollison [19], Thieme [25], Thieme and Zhao |26] for some examples. 

However, the nonlinearities in equations/systems often give expression to the inter- or 
intra-specific actions in population dynamics. Intuitively, the effect of nonlinearities should 
be reflected by many dynamical properties including the asymptotic speeds of spreading. 
Namely, the linear determinacy of minimal wave speed of traveling wavefronts and asymp- 
totic speeds of spreading cannot be true for all nonlinear models. For autonomous scalar 
equations, a famous counter example of linear determinacy is 



(0,0), (1,0), (0,1) 



and K = (k\, k 2 ) defined by 
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where v > — 1 is a constant that does not appear in the following linearized system 



and we refer to Hadeler and Rothe [9] for precise results on its asymptotic speed of spread- 
ing, which is not linearly determinate for v > 2. Moreover, some results on asymptotic 
spreading have also been obtained for coupled diffusion systems with multi equilibria, 
which formulates the role of inter-specific coupled nonlinearities, see Lin et al. |15j and 
Weinberger et al. [28] for two examples of integral-difference equations, Lin |13| for a 
predator-prey reaction-diffusion system. 

For reader's convenience, we first give the following definition. 

Definition 1.1 Assume that u(t,x) is a nonnegative function for x E K, t > 0. Then c* 
is called the asymptotic speed of spreading of u(t, x) if 

a) lim^oo sup| i .| >(Ci+e)t u(t,x) = for any given e > 0; 

b) liminf^oo inf^^.^ u(t,x) > for any given e G (0,c*). 

Clearly, the asymptotic speed of spreading states the observed phenomena if an observer 
were to move to the right or left at a fixed speed [27]. Biologically, it also describes the 
speed at which the geographic range of the new population expands [10]. Therefore, it 
becomes a very important index formulating the spatial propagation of ecological com- 
munities. At the same time, it is possible that the asymptotic speed of spreading of a 
nonnegative function is not a positive constant in the above limit sense, see Berestycki 
et al. [3] for some examples. When the asymptotic speed of spreading is not a constant, 
its lower bounds and upper bounds in [21 Section 1.8] are still useful because these can 
describe and estimate the success of biological invasions. 

If an irreducible cooperative system has just two equilibria in the interesting inter- 
val, it is very likely that all the unknown functions have the same asymptotic speed of 
spreading coincided with the linear determinacy, see some results by Liang and Zhao [121 . 
Lui [ISl [17]. In particular, when (jl.ip is concerned, Li et al. [HI Example 4.1] studied 
the propagation modes when one species is the aboriginal and the other is the invader, 
namely, the interesting interval is 



on which the system has no other equilibria, and can also be studied by [T^USIIIT] . 

In this paper, we consider the asymptotic spreading of (jl.ip when both species are 
invaders, namely, (0, 0) will be the invadable equilibrium and the interesting interval will 




Au(t, x) + u(t, x), 



at 



[Mi] x [0,fc 2 ] or [0,A;i] x [l,fc 2 ] 
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be [0, ki] x [O,/^]) on which (jl.ip has four equilibria such that we cannot use the theory 
of [T2j HUE]. To obtain some estimates on asymptotic spreading, the abstract results 
developed by Berestycki et al. [2] will be applied, and the lower bounds of asymptotic 
speeds of spreading will be estimated. More precisely, we first give some properties of u\, 
then we regard the second equation of (jl.ip as a nonautonomous equation and establish 
some conclusions by [2]. Our results imply that: (1) The nonlinearities described the 
inter-specific actions may play an important role in asymptotic spreading such that the 
asymptotic spreading of one species is faster than the case that the inter-specific actions 
disappear; (2) It is necessary to use different indices to formulate the asymptotic spreading 
of each unknown functions if the system has multi equilibria. Moreover, our results answer 
the nonexistence of traveling wave solutions of (jl.ip . which also develops the theory of 
traveling wave solutions in Lin et al. |14j . 

In Section 2, we shall give some preliminaries, including a classical conclusion of Fisher 
equation and an important result established by Berestycki et al. [2]. Then we shall show 
some estimates on the asymptotic spreading of (jl.ip if both species are invaders, which 
are also applied to the study of the corresponding traveling wave solutions. In the last 
section, further discussion is provided to illustrate our conclusions. 

2 Preliminaries 

We first present some results of the following Fisher's equation 



in which all the parameters are positive and z(x) > is a uniform continuous and bounded 
function. Due to the theory of asymptotic spreading established by Aronson and Wein- 
berger p], we have the following result. 




dAz(t, x) + rz(t, x) [1 - z(t, x)/K] 
z(x), 



(2.1) 



Lemma 2.1 Assume that z(t,x) is defined by (|2.ip and e £ (0, 2\fd~r) holds. Then 



lim 



inf 

\x\<(2Vdr-e)t 



z(t,x) = K. 



Moreover, if z{x) admits compact support, then 



lim 



sup 

\x\>{2Vdr+e)t 



z(t,x) = 0. 



For (|2.1j) . the following comparison principle is also true (see Ye and Li [29j). 
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Lemma 2.2 Assume that z(t, x) > 0, x G R, t > 0, satisfies 

^f 1 > (<)dAz(t, x) + rz(t, x) [1 - z(t, x) /K] , 

z(0,x) > (<)z(x). 

Thenz(t,x) > (<)z(t,x), where z(t,x) is defined by (|2.ip . 

For the system (jl.ip . we also give the following comparison principle (one also refers 
to Pao [21], Smoller [24] . Ye and Li |29] for more details). 

Lemma 2.3 Let (ui(t,x),u 2 (t,x)) be defined by 

dui ^ = diAui(t,x) + riui(t,x) [1 - ui(t,x) + b x u 2 (t,x)] , 
du2 ^ = d 2 Aii 2 (£, x) + r 2 n 2 (t, x) [1 - u 2 (t, x) + 6 2 iii(i, x)] , 
«i(0,x) = ui(x),u 2 (0, x) = u 2 (x), 

where U\{x) > 0, n 2 (x) > are uniformly continuous and bounded. If (zi(t, x), z 2 (t, x)) > 
(0, 0) is uniformly continuous and bounded for (i, x) G (0, +oo) x R and satisfies 

> (<)d 1 Az 1 (t,x) + r 1 z 1 (t,x) [1 - z x {t, x) + b x z 2 (t, x)] , 

> (<)d 2 Az 2 (t,x) + r 2 z 2 (t,x) [l-z 2 (t,x) 
^i(0,x) > {<) Ul {x),z 2 (0,x) > (<)n 2 (x). 

T/ien (z!(t,x),z 2 (t,x)) > (<) (ni(t, x), n 2 (t, x)). 

Now, we consider a nonautonomous equation as follows 

dAu(t, x) + f(t, x, n), 
u{0, x) = n(x), 

in which / : RxRxR + — > R is assumed to be of class C" 5 / 2 '* 5 in (t, x), locally in u, for a given 
5 € (0, 1). Moreover, / is also locally Lipschitz continuous in u and of class C 1 in u G [0, /3] 
with /? > uniformly with respect to (i, x) G R x R, it is also supposed that f(t, x, 0) = 0. 
Since (|2.2p cannot generate a semiflow, the study of its asymptotic spreading is very hard. 
To formulate its asymptotic spreading, we first present some important definitions and 
results given by Berestycki et al. [21 Section 1.5]. 

Definition 2.4 We say that complete spreading occurs for a solution u(t, x) of (|2.2p if 
there is a function t — > r(t) > such that r{t) — > oo as t — > oo and the family (B r ^) t yQ 
is a family of propagation sets for u, that is 



b 2 z 1 (t,x)\ , 



du(t,x) 

at 



(2.2) 



lim inf 

t— >oo 



inf x) > > 0, 
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where B r = {x € R : \x\ < r}. 

This definition, in fact, gives a description of the success of spatial spreading/invasion, 
which is similar to the second item of Definition 11.11 Since there are only two directions 
in R, we also show a specific case of Berestycki et al. [21 Definition 4] as follows. 

Definition 2.5 We say that a family (r(i))t>o of nonnegative real numbers is a fam- 
ily of asymptotic spreading radii for a solution u(t, x) of (|2.2p if the family of segments 
([— r (£), r(t)])t>o is a family of propagation sets for u(t,x), that is 

lim inf < inf u(t,±s) > > 0. 
t^oo Lss[0,r(t)] J 

Definition 2.6 We say that a family (r(t))t>o is a family of admissible radii if (r(t))t>o £ 
C 1 + S / 2 (R+,R+) and sup t > \r'(t)\ < oo. 

Remark 2.7 If lim inf t->oo r (t)/t exists and is positive, then it is a lower bounds of asymp- 
totic speed of spreading. Such a definition is still useful because it can describe the 
phenomena of successful invasion, even if lim^oo r(t)/t does not exist, see Berestycki et 
al. [21 Section 1.8] for the upper bounds of asymptotic speed of spreading. 

For (j) £ C X ' 2 (R x M), define 




-dA<t>-f u {t,x,Q)<l>. 



Considering the generalized principal eigenvalue problem formulated by 

A; = inf{A € R,3(p £ C 1>2 (R x R) Q W 1,OQ (R x R), inf (f> >0,L<j)< A0}, 
then A' x < implies that the equilibrium is unstable and the following conclusion holds. 



Lemma 2.8 ( |2j) Assume that A' x < and there exists r(t) of admissible radii such that 

lim inf u(t, ±r(t)) > 0. 



t— ¥00 

Then 



t— >oo 



lim inf <^ inf u{t,x) } > 0. (2.3) 



J*|<r(t) 

Lemma 2.9 ( [2j) Letu(t,x) be the solution of the Cauchy problem (]2.2p associated with 
an initial datum u{x) > 0. Assume that A' x < holds and there exists r(t) of admissible 
radii such that 

lim inf 1 lim inf { inf (Mf'(t, x ± r(t), 0) - (r'(t)) 2 ) 1 1 > 0. (2.4) 

R-^oo [ t^tco Y\x\<R J J 

Then (|2.3|) holds for u(t,x). 
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3 Main Results 

In this section, we first prove the following result on asymptotic spreading. 

Theorem 3.1 Let (ui(t,x),u 2 (t,x)) be defined by 

&Ul ^ = d%Aui(t, x) + ri«i(t, x) [1 - ui(t, x) + b%u 2 (t, x)] , 
9u f t ,x) = d 2 Au 2 (t, x) + r 2 u 2 (t, x) [1 - u 2 (t, x) + b 2Ul (t, x)\ , (3.1) 
^ui(0,x) = </> 1 (x),u 2 (0,x) = (fo(x), 

in which <fii(x) > 0,(j) 2 (x) > are uniformly continuous and bounded for Suppose 
that d\r\ > d 2 r 2 holds. Then 

lim inf u 2 (t,x) = lim sup u 2 (t,x) = k 2 (3.2) 



and 



lim inf u\(t,x) = lim sup u\(t,x) = k\ (3-3) 

t^co\ x \ <c t ^°°\x\<ct 



for any c < c* = min{2 v / c£iri, 2^ d 2 r 2 {\ + b 2 )}. 

For the main condition of the theorem, we give the following remark. 

Remark 3.2 By [1J, d\r\ > d 2 r 2 implies that u\ has stronger spreading ability than that 
of u 2 if the inter-specific actions disappear in (|3.ip (namely, b\ = b 2 = in 



Before verifying Theorem 13.11 we first prove several lemmas, through which the con- 
ditions of Theorem 13.11 will be imposed. 

Lemma 3.3 For all x € R, t > 0, the Cauchy problem (|3.1j) admits a unique solution 
(ui(t,x),u 2 (t,x)) such that 

(0,0) < ( Ul (t,x),u 2 (t,x)) < {Ex,E 2 ), 

in which 

Ei = max {sup 01 (x), ki, ^-supcp 2 (x)\ ,E 2 = max {sup <j) 2 {x), k 2 , ^ sup 4>\ (x) 1 . 

liceK k 2 x£ k J IxER k\ x gK J 

Proof. We prove the lemma by comparison principle. Clearly, 

n#i(l - E x + biE 2 ) < 0,r 2 E 2 (l -E 2 + b 2 E x ) < 0. 
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Then (Ei,E 2 ) is an upper solution while (0,0) is a lower solution of (|3.ip for (t,x) G 
(0, +00) x R. Therefore, Lemma 12.31 indicates that 

(0,0) < { Ul (t,x),u 2 (t,x)) < (E 1 ,E 2 ), (t,x) e (0,+oo) x R. 

The strict inequalities are evident by the following two facts: 

(1) The heat operator has property of infinite propagation speed; 

(2) (pi(x),(j)2(x) admit nonempty supports. 

The proof is complete. □ 

Lemma 3.4 Define c\ = 2\Jd\r\. Suppose that u\{t,x) is defined by (|3.ip . Then 



liminf inf u\(t,x) > 1 

t-5>oo \ x \<ct 



for any c < c\. 



Proof. By Lemma 13.31 we see that 
dui(t, x 



dt 



> diAui(t, x) + rmi(t, x) [1 - m(t, x)] , x G R, t > 0. 



Then the result is evident by Lemmas 12.11 and 12.21 The proof is complete. 
Let /3 > be a constant such that 

Pm + ri«i [l — ui + biu 2 ] , f3u 2 + r 2 u 2 [1 — u 2 + b 2 u{\ 

are monotone increasing if 

(0,0) < ( Ul ,u 2 ) < (E U E 2 ). 
For t > 0, define T(t) = (Ti(t), T 2 {t)) as follows 



Ti(t)ui{x) 



0-y) 2 

it poo — — / \ 

VI -=j_ oo e « lt Ul{y)dyi 



-Pi 



[ T *®M*) = IZo e-Vu 2 (y)dy. 

Moreover, for t > 0, s > 0, we still denote 



T l ( S > *) = ^ Ho e " 4dlt 2 "1 ( S > V) d V> 

T 2 (t)u 2 (s,x) = -0^J^ o e-^u 2 (s,y)dy. 



(3.4) 



□ 
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Let X be defined as follows 

X = {u : u is a bounded and uniformly continuous function from WL to M 2 }, 

which is a Banach space equipped with the supremum norm. Then T(t) : X — > X is an 
analytic semigroup (see [B]). Denote 

X + = {u : u 6 X,u > 0}. 

Then T(i) : X + — )■ X + is a positive semigroup. Using the standard theory of semigroup 
(see |22|). we have the following conclusion. 



Lemma 3.5 The unique solution of ()3. 1|) can also be formulated by 

Is. 

(3.5) 



ui(t,x) = Ti(t)(f)i(x) + f Q Ti(t - s)[Fi(ui,u 2 )]{s,x)ds 1 



u 2 (t,x) = T 2 (t)(p 2 {x) + foT 2 (t - s){F 2 (u 1 ,u 2 )](s,x)ds, 

in which Fi(ui,u 2 ) = j3u\ + riui[l — u\ + hu 2 ], F 2 (u\,u 2 ) = (3u 2 + r 2 u 2 [l - u 2 + b 2 u{\. 

By above lemmas, we give the proof of Theorem 13.11 as follows. 
Proof. Let c < c* be fixed. By (|3.4p . we can choose e > satisfying the following facts. 



(A) There exists T > such that 



inf m(t, x) > 1 - e for all t > T. (3.6) 

4|x|<(c+3c*)t 



(B) 4^2^(1 + 62(1-6)) > c + c* > 2c. 

Define r(t) = (c + c*)t/2. Then r(t) G C°°(IR + ,1R + ) such that Definition [221 is true. 
Moreover, lim^oo r(t) = 00 also implies that Definition [23] holds and a complete spreading 
of u\ has been proved. 

Denote 

du 2 (t,x) 



Ot 



d 2 Au 2 (t, x) + r 2 u 2 (t, x) [1 - u 2 (t, x) + b 2 ui(t, x)] 
: d 2 Au 2 (t,x) + f(t,x,u 2 ), (3.7) 



in which the definition of / is clear. To apply Lemma 12.91 we encounter some difficulties 
since f(t, x, u 2 ) has no definition if t < 0. So we define / such that 



f(t,x,u 2 ) 



J(t,x,u 2 ),t > 1, 
g{t,x,u 2 ),t e [0,1], 
r 2 u 2 (l - u 2 ),t < 0, 



10 



G. Lin 



in which g(t,x,U2) = u 2 g%(t,x,U2) with 

r 2 (l - u 2 ) < gi(t, x, u 2 ) < r 2 (l - u 2 + b 2 ui),t G [0, l],x G R 

such that / satisfies the smooth condition of (|2.2|) . Since [0, 1] is a bounded interval, the 
existence of / or g is clear. Consider the following initial value problem 

^El = d 2 Az 2 (t,x)+J(t,x,z 2 ) 
^ z 2 (0,x) = z{x) > 0, 

and 

= d 2 Az(t, x) + f(t, x, z), 
z(0,x) = z(x) > 0. 
Then the comparison principle implies that 

z 2 (t,x) > z(t,x), (t,x) £ (0,oo) X 



(3.8) 



Thus, it suffices to study 

du 2 (t,x) 



dt 



d 2 Au 2 (t, x) + f(t, x, u 2 ). 



Evidently, for any (i, i) 6 I x I, we have 

r 2 E 2 > fi 2 {t,x,0)>r 2 . 

For the Fisher equation 

du 2 (t, x) 



dt 



d 2 Au 2 (t, x) + r 2 u 2 (t, x) [1 — u 2 (t, x)\ , 



(3.9) 



we see that the corresponding X[ < by Berestycki et al. (2j Section 1.5]. Then (]3.9 
implies that the corresponding X[ of (|3.7f) is also negative. 

Therefore, we just need to verify that (12. 4p is true. For any R > 0, we see that 



hminf inf [4< 2 (t, x ± r(i), 0) - (r'(t)) 2 ] 

t^oo {\x\<R L 



liminf < inf 

t^oo \x\<R 



4df'(t,x±r(t),0) 



c + c* 



liminf< inf [4dr 2 (l + b 2 u\(t, x ± r(t)))\ 
t^-oo \x\<B. 



C + C* 
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By the item (B), it is clear that (|3.6p holds if t > is large enough. Therefore, 
liminf { inf [Adr 2 (l + b 2 u 1 (t,x ±r(t)))} - ( \ 

i->oo l|x|<R \ 2 J \ 

> 4dr 2 (l + 6 2 (l- e ))-(^) 2 

> 

since R < —^-^-t — > oo as t — > oo. 

Note that 4dr 2 (l + 62(1 — e)) — ls independent of R > 0, we also obtain that 

liminf jliminf j inf \Uf' (t,x ± r(t), 0) - (r'(t)) 2 ] \ \ > 0. 

_R-s>oo [ t->oo (ja;|<.R J J 

By Lemmas ESES and (pT8|) . if c = c, then 

liminf inf it2(i, x) > 0. 
t^oo |a|<ct 

Due to the arbitrary of c, what we have done implies that 

liminf inf u 2 (t,x) > 0, liminf inf u\(t,x) > 

for any c < c*. 

It suffices to verify that ()3. 2f) and (13. 3ft are also true for any fixed c < c* . Let 2c\ = 
c + c* , c\ > C2 > • • • > c n > c n+ \ > • • • , linin^oo c n = c, define positive constants 

liminfi^ 00 inf| :! ,|< Cnt -ui(>,x) = uf, lim'mf t - ¥oo m.i\ x \ <Cnt u 2 (t,x) = u 2 , 
i limsup^ 0O sup| iB | <Cnt «i(t,x) limsup t _ >00 sup| a .| <Cnt w 2 (i,a;) = u 2 l - 

and 

liminf t _ ) . 00 inf| a .| <3 ui(t,x) = u x , liminf^oo inf^o n 2 (t, x) =u 2 , 
^limsup t _ >00 sup| a! | <a ui(t,x) =«i, limsup t _>. 0O sup| a .| <Et M 2 (t,a;) = u 2 . 

Clearly, these positive constants are well defined and satisfy 

(LI) ;u™,u 2 are nondecreasing and u% < u.i,u 2 < u 2 for all n > 0; 

(L2) u™,«2 are nonincr easing and > %, > u 2 for all n > 0; 

(L3) lim^-Kxju" and lim„_ ) . 0Q u^ 1 exist for i = 1,2; 

(L4) lim^oo < Ui <Ui < lim n _>oo , « = 1, 2. 
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For each n > 1, t — > oo implies that (c n +i — c n )t — > oo. Using the positivity of the 
semigroup of T{t) and the dominated convergence theorem for t — > oo in (j3.5j) . we see that 

_ n+ l . ffff +r 1 rff[l-rff + feii%] 

by the monotonicity of F\. Letting n — > oo, we further obtain that 

1 - lim ul + 6i lim u% > 0. 



In a similar way, we have 



1 - lim u™ + bi lim jxg < 0, 

n— >oo n— >oo 

1 - lim + 6 2 lim u% < 0, 

n— >-oo n-^-oo 

1 - lim «2 + b 2 lim H™ > 0, 



and 

lim = lim H™ = fei, lim U2 = lim H?? = /C2. 

n— >oo n— >oo n— >oo n— >oo 

From (L4), we obtain 

Since c is arbitrary, we complete the proof. □ 
We now present three remarks to further illustrate our conclusion. 

Remark 3.6 If d\ = <i 2 ,ri = then Lin et al. [14] implies that (II. lh has a traveling 
wave solution connecting (0,0) with (fci , ^2) for any wave speed which is larger than 
1\fd~\r\ = 2y 'd2T2- Therefore, if < 4>i{x) < &i,0 < <p2{x) < ^.2 admit compact supports, 
then the standard comparison principle states that the asymptotic speeds of spreading of 



two invasion species are not larger than 2y/d\r\ (see the subsequent Propositions 13.91 and 
13. lQf) . By Lemma 13.41 the asymptotic speeds of spreading of both invasion species are 
2^d\r\ = 2^d 2 r 2 . 

Remark 3.7 If d\T\ > c^^/^ = (^2(1 + &2&1) holds and < <j)\{x) < k±,0 < </> 2 (x) < k,2 
admit compact supports, then < Ui(t,x) < ki and 

8U2(t, X) , . . 1 \ r j / \1 

— < d 2 An 2 (t,x) + r 2 U2{t,x) [k 2 - u 2 {t,x)\ . 

at 

Namely, « 2 is a lower solution of the following Cauchy problem 

dw2 ^ x ) = (J 2 Aw 2 (t, x) + r 2 u> 2 (i, x) [k 2 - W2(t, x)] , 
w 2 (0,x) = 4> 2 (x). 
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Then the comparison principle (Lemma 12 .2|) indicates that v,2(t,x) < w 2 (t,x), and the 
upper bounds of asymptotic speed of spreading of U2(t, x) is not larger than 2 \/ 0(2^2^2 by 
Lemma 12.11 Recalling Lemma 13.41 the lower bounds of asymptotic speed of spreading of 
u\(t,x) is larger than 2\fd\r\ such that two species have two distinct asymptotic speeds 
of spreading even if both of them are constants. 



Remark 3.8 If 2^d 2 r 2 < 2\[d^F\ < 2y/d 2 r 2 {l + b 2 ) with d\ = d 2 and < <j>i(x) < 
k%, < (p2(x) < k 2 admit compact supports, then Lin et al. [H] implies that the asymptotic 
speeds of spreading of both invasion species are less than 2^/d\r% (see Propositions 13.91 and 
13. lOf) . and Theorem 13. II indicates that the asymptotic speeds of spreading of both species 



are 2y/d\r\ such that the invasion of 1x2 is fastened by u\. 

Before ending this section, we also apply our main result to the study of traveling wave 
solutions of ([Lip . 

Proposition 3.9 If (jl.ip has a traveling wave solution (ux(t,x),u 2 (t,x)) = (ipi(x + 
ct),tjj2(x + ct)) connecting (0,0) with (kx,k 2 )- Then the asymptotic speeds of spreading 
of u%(t, x), u 2 (t 1 x) are not larger than c if < cfri(x) < k±,0 < 4> 2 {x) < &2 admit compact 
supports. 

Proof. In the lemma, a traveling wave solution (ui(t,x),U2(t,x)) = (ipi(x+ct),i(j2{x+ct)) 
connecting (0,0) with (^1,^2) is formulated by 

lim (Vi(s),V 2 (s)) = (0,0), lim (Vi(s),V 2 (s)) = (h,k 2 ), 

s— >— OO s— ^00 

where (ipi,ip2) is the wave profile and c is the wave speed. 

For any p G R, a traveling wave solution (ipi (x + ct + p), ip2 (x + ct + p)) is also an entire 
solution (defined for all t & M.) of the following Cauchy problem 

d"^'^ = d\Aui(t,x) + riui(t, x) [1 - ui(t,x) + b\u 2 {t,x)\ , 
9u2 tt' x ) = d 2 Au2(t,x) + r 2 U2(t,x) [1 - u 2 (t,x) + 62^1 (t,x)] , 
ut(0,x) = V>i(x + p),u 2 (0,x) = ip 2 (x + p). 
Letting p large enough, then 

ipi(x + p) > (f>i(x),ip 2 (x + p) > (j) 2 (x) 

since < (f>i(x) < k%, < 4>2{x) < k2 have compact supports. Now (-01 (x + ct + p) , ip2 (x + 
ct + p)) becomes an upper solution of (13. ip . Then the comparison principle implies that 

ipi(x + ct + p)> ui(t,x),i/i 2 (x + ct + p) > u 2 (t,x) 

and the result is clear. □ 
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Proposition 3.10 Under the assumptions of Remark \3.6\ or Remark \3.8l (jl.ip has a trav- 
eling wave solution connecting (0,0) with (ki,k 2 ) if c > m&x{2\/ d\r±, 2\/d 2 r 2 }. Moreover, 
if di > rf2> r i > r 2> i/ien (jl.ip /ias a traveling wave solution connecting (0,0) with (kx,k 2 ) 
if c > 2\fd\r\. If c < 2y/d\r\, then (jl.ip /10s no£ a traveling wave solution connecting 
(0,0) with (fa,k 2 ). 

Proof. If c > max{2\/G£iri, 2\/^2^2}j then define < 7^1 < 7i2 by 

^7ji - Ha +ri = daf 2 - C7i 2 + n = for i = 1, 2. 

By Lin et al. [HI Theorem 5.11], the lemma is true if (711,712) 0(7217 722) is nonempty, 
which is evident if Remark 13.61 or Remark 13.81 holds or d\ > d 2 and 77 > r 2 are true. 

The nonexistence of traveling wave solutions is clear by Lemmas 12.11 and 13. 4[ and we 
omit the proof here. The proof is complete. □ 

4 Discussion 

In ecological systems, the cooperatitive/symbiotic/mutualistic communities are very 
universal. For example, the role that insects, in particular bees, have in the fecundation 
of flowers, see Boucher [5]. Furthermore, Malchow et al. |18j Section 4.3.2] also introduced 
many examples. In population dynamics, the behavior of many cooperative kinetic systems 
is very simple: If a cooperative system admits only one positive equilibrium, then the 
equilibrium is asymptotic stable and the others are unstable. These mathematical results 
are very easy and can be found in many textbooks, we also refer to Malchow et al. [18], 
Murray [20] . In particular, if 6162 < 1 in (jl.ip . then (£7, k 2 ) is stable and the phase plane 
of the corresponding kinetic system is very clear, see Murray |20j pp. 101]. Biologically, 
(ki,k 2 ) > (1,1) implies that each species has increased its steady state population from 
its maximum value in isolation [20], which is achieved by inter-specific cooperation. 

However, when the spatial-temporal structure is involved in cooperative systems, e.g., 
the spatial dispersal of plant and seeds (see Murray [2D1 Section 3.6]), its dynamical prop- 
erties may be very complex since the process often involves the far-from-equilibrium dy- 
namics. By Liang and Zhao [12j . Lui |16jll7j. if an irreducible cooperative system admits 
two steady states, it is very likely that different unknown functions have the same asymp- 
totic speed of spreading. However, Remarks I3.6H3.8I show that the complex propagation 
modes of evolutionary systems with multi equilibria since it is necessary to formulate the 
asymptotic spreading of different unknown functions by different indices. Note that the 
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number of steady states is determined by the nonlinearities, this certainly indicates the 
complex arising from the nonlinearities. 

We now consider the linear determinacy problem. Because we consider the spatial 
invasion of two species, then one interesting equilibrium is (0, 0) that is invadable. If the 
asymptotic speeds of spreading are linearly determinate, then the asymptotic speeds of 
spreading will be fully determined by d\, n, d%, T2, which is impossible by Remarks l3.6H3.8l 
Therefore, our results show the effect of inter-specific cooperation from the following two 
factors: (1) asymptotic speed of spreading or its lower bounds of U2 since c* > 2\fd^X2 
in Theorem 13 .1\ (2) eventual population densities on the coexistence domain because of 

In this paper, utilizing the theory established by Berestycki [2], we obtain some esti- 
mates of the asymptotic speeds of spreading, which partly shows the role of nonlinearity. 
Unfortunately, only the lower bounds and upper bounds of asymptotic speeds of spreading 
are obtained, precise results need further investigation. 
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